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ABSTRACT 

It is usually believed that the Cowling approximation can give satisfactory solutions 
if the stars have soft equations of state and/or if the strongly general relativistic 
stars are treated in the case of rigid rotation. Since, however, there have been no 
systematic studies about the accuracy of the Cowling approximation for differentially 
rotating compressible stars, we investigate eigenfrequencies and eigenf unctions of the 
oscillation modes in rapidly and differentially rotating compressible stars by employing 
the exact method including full-perturbations and the Cowling approximation. 

We have found that the Cowling approximation for f-mode oscillations is not a 
good approximation in rapidly and differentially rotating stars, although rapid rotation 
makes this approximation better for rigidly rotating stars. This result suggests that we 
must be careful when we apply the Cowling approximation to differentially rotating 
stars even in the framework of general relativity. On the other hand, the approximation 
will work well for r-modes even if the star is rotating differentially. Therefore, the 
Cowling approximation can be used as a strong tool for the investigation of r-mode 
oscillations in the general relativistic framework that it is difficult to compute including 
the perturbations of gravity. 

Key words: gravitation - instabilities - stars: oscillation - stars: rotation. 



1 INTRODUCTION 

Concerning the oscillations of spherical stars. Cowling 
(1941) started to use an approximation that neglect the per- 
turbed gravitational potential to rewrite the system of equa- 
tions into easier form. By using this approximation (later it 
is called the Cowling approximation), he found that there 
are three types of oscillations; the high frequency acoustic 
mode (p-mode), the low frequency gravity mode (g-mode), 
and the intermediate fundamental mode (f-mode). 

In the last two decades, the Cowling approximation has 
been applied to general relativistic compact stars. In spheri- 
cal stars, McDermott, Van Horn, and SchoU (1983) proposed 
the scheme in which all components of the metric perturba- 
tions are neglected. On the other hand, Finn (1988) pro- 
posed the formulation which neglects most components of 
the metric perturbations except for the (r, t) component. 
The accuracies of these two approaches were compared by 
Lindblom and Splinter (1990) for mainly p-modes, and it 
was found that generally the approximation neglecting all 
components of the metric perturbations is better than the 
other one, while there is a possibility that the latter be- 
comes better when the g-modes are considered. Hereafter, 



"the Cowling approximation" means the one that neglects 
all components of the metric perturbations. 

It is also shown that the Cowling approximation is re- 
liable in slowly rotating general relativistic stars (Kojima 
1997). Furthermore, Ipser and Lindblom (1992) proposed 
the formulation of applying the Cowling approximation to 
rapidly rotating relativistic stars. The accuracy of the Cowl- 
ing approximation in rotating stars was studied by Yoshida 
and Kojima (1997). They compared the two sets of solu- 
tions obtained by using the procedure produced by Ipser 
and Lindblom (1992) and by solving the exact equations 
for the f- and p-mode oscillations of nearly spherical stars 
with slow rotation. According to their results, solutions of 
the Cowling approximation agree well with the exact so- 
lutions for slowly rotating and highly relativistic models. 
Recently, Yoshida and coworkers have used the Cowling ap- 
proximation to find neutral points of the f-mode oscillations 
of rapidly rotating polytropes (Yoshida and Eriguchi 1997, 
1999, 2001; Yoshida, RezzoUa, Karino, and Eriguchi 2002). 
Also they checked the accuracy of the approximation in rel- 
ativistic stars and found the conditions that the approxima- 
tion becomes accurate (Yoshida and Eriguchi 1997, 1999). 

Our knowledges about the Cowling approximation to 
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the oscillations of rotating stars thus far can be summarized 
as follows: 

The Cowling approximation becomes more accurate 

(1) when the compressibility of the star is sufficiently high, 

(2) whcn the gravity of the star is sufficiently strong. 

It may be natural that the Cowling approximation becomes 
more accurate if the effect of the gravity does not work effec- 
tively to the fluid mass. For the simple oscillations such as 
f-mode, the amplitude of the mode becomes a simple mono- 
tonically increasing function of a stellar radius, and takes its 
maximum at the surface. Here, when we consider the density 
distribution of a compressible star, the density of surface re- 
gion becomes relatively low. In other words, there are only 
small masses in the region where the mode amplitude gets 
large, in the centrally condensed star. Also, effectively, the 
strong gravitation in general relativity (GR) can have the 
same effect. Consequently, the effects of compressibility and 
strong gravity make the accuracy of the Cowling approxi- 
mation better, 
and 

(3) when the rotation of the star is sufficiently rapid, when 
the star rotates rigidly. 

Here, notice that the effect of differential rotation on the 
approximation has not studied, therefore we will focus on 
it. 



From the astrophysical points of view, we should take 
additional effects into account. For example, some stars 
are rotating differentially. In particular, newly born neu- 
tron stars may rotate differentially (see e.g. Liu and Lind- 
blom 2001, Shibata and Uryu 2000, Dimmelmeicr, Font and 
Miiller 2002). Young neutron stars have been considered 
to be astrophysically interesting objects because they can 
be candidates of sources of gravitational waves due to the 
rotational instability (Chandrasekhar 1970, Friedman and 
Schutz 1978ab). However, if stars rotate differentially, the 
eigenfunctions of oscillations as well as the equilibrium con- 
figurations will be deformed considerably (see e.g. Eriguchi 
and Miiller (1985) for Newtonian equilibrium stars; Ko- 
matsu, Eriguchi, and Hachisu (1988) for relativistic equilib- 
rium stars; Karino, Yoshida, and Eriguchi (2001) for eigen- 
functions of oscillations of Newtonian staxs.) If one would 
like to apply the Cowling approximation to such differen- 
tially rotating stars, it is necessary to estimate the applica- 
bility of the Cowling approximation to them. In this paper 
we will study the accuracy of the Cowling approximation 
in rotating Newtonian stars by comparing the results ob- 
tained by the Newtonian Cowling approximation with those 
by solving the exact equations for the perturbed equations. 



The plan of this paper is as follows. In the next sec- 
tion, we will briefly summarize the scheme to solve the lin- 
earized perturbed equations for differentially rotating stars 
with and without the gravitational perturbation. Our nu- 
merical results will be presented in Section 3. In Section 
4, the accuracy of the Cowling approximation will be dis- 
cussed. In the last section, we will summarize our present 
study and make a brief comment on the applicability of the 
Cowling approximation to GR. Throughout this paper, the 
ordinary spherical polar coordinate system (r, 9, if) is used. 



2 METHOD 

The motivation of this paper is to estimate the accuracy 
of the Cowling approximation. For this purpose, we will 
solve the linearized equations of perturbed fluid in New- 
tonian gravity and numerically obtain the eigenvalues and 
eigenfunctions of oscillation modes. Then we will compare 
them with the results computed by neglecting the gravity 
perturbation, Scf), from the basic equations (the Cowling ap- 
proximation) . The method of solving the full system of the 
basic equations for the oscillations of rotating stars is shown 
in Karino, Yoshida, Yoshida and Eriguchi (2000). Hence, in 
this section, we will briefly explain the solving method. 

The equilibrium states of axisymmetric polytropes can 
be obtained by the method proposed by Eriguchi and Miiller 
(1985). The polytropic relation is expressed as 

p = Kp'+^, (1) 

where p, p, K and A'^ are the pressure, the density, the poly- 
tropic constant and the polytropic index, respectively. In 

this paper, we compute polytropes with A'^ = 0.5, 1.0, and 
1.5. Such relatively stiff polytrope can be used as the models 
of neutron stars. 

We will use the following rotation law (j-constant rota- 
tion law): 

o- (2) 
(rsine)2+^2' 

where SI and Oo arc the angular velocity and the angular 
velocity on the rotation axis. The parameter A represents 
a degree of differential rotation. A small value of A corre- 
sponds to a strongly differential rotation, while in the limit 
of ^ ^ oo the rotation law becomes that of rigid rotation. 

The eigenfrequencies and eigenfunctions of the modes 
can be obtained by solving the perturbed fluid equations. 
Here, we assume small perturbations and linearize the sys- 
tem of equations. The linearized quantities are expanded as 
follows: 

Sfir,e,<p,t) = ^exp{i{at - m<p))fm{r,e), (3) 

m 

where / and 5f are a certain physical quantity and its Eule- 
rian perturbation, respectively. In this expansion, cr denotes 
the eigenvalue of the oscillation mode and m is the azimuthal 
wave number. Here, we only consider the rn — 2 i- and r- 
mode for two reasons: (1) from the Eistrophysical point of 
views, the m = 2 mode is the most interesting because it 
has the fastest growth rate in rapidly rotating stars, and (2) 
it has been said that the accuracy of the Cowling approxi- 
mation becomes the worst for the m = 2 mode. Actually, the 
eigenfrequencies of radial and quasi-radial-modes have been 
studied by nonlinear computations and it has been shown 
that the frequency obtained by the Cowling approximation 
is considerably different from the one given by full GR calcu- 
lations even in the case with rapid rotations (Font, Goodale, 
Iyer, Miller, Rezzolla, Seidel, Stergioulas, Suen and Tobias 
2002). 

On the other hand, in the Cowling approximation, 
the Eulerian perturbation of the gravitational potential is 
dropped from the full system and solved the equations with 
the remaining terms. 
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3 NUMERICAL RESULTS 

The eigenvalues of the m = 2 f-modes for the N = 1.0 poly- 
trope are plotted in Fig. 1 as a function of T/|W^|, where T 
is the rotational kinetic energy and W is the gravitational 
potential energy of the star. Three curves of the different 
types correspond to the different values of the parameter 
A; A = 0.5, 1.0, and oo (= rigid rotation). For each A, two 
curves are plotted: the upper curves correspond to the re- 
sults of the full computations and the lower to the approx- 
imated solutions. In the case of ^ = 1.0 and 1.5 including 
the perturbation of gravity, when the rotation of the star get 
rapid, the so-called neutral point where a = appears. This 
neutral point is corresponding to the critical point of secu- 
lar instability due to gravitational radiation reaction (Chan- 
drasekhar 1970, Friedman & Schutz 1978ab). On the other 
hand, the curves given by the Cowling approximation have 
smaller, finite values. This means that we can not estimate 
the critical points of secular instability when we apply the 
Cowling approximation. 

The differences of the eigenfrequencies between the re- 
sults of the full perturbed systems and the Cowling approx- 
imation are shown in Fig. 2. From this figure, it is clear that 
the accuracy of the Cowling approximation becomes better 
as the stellar rotation gets faster for rigidly rotating poly- 
tropes. For differentially rotating polytropes, however, the 
accuracy of the approximation becomes worse and worse as 
the stellar rotation becomes faster. 

The comparison between the approximated results and 
full perturbed results have been done for rigidly rotating 
polytrope by Jones, Andersson and Stergioulas (2002). They 
have carried out time evolving computations of the lin- 
earized perturbed equations of rapidly and rigidly rotating 
Newtonian polytrope with the Cowling approximation, and 
compared their results about f- and r-modes with the previ- 
ous results. They have checked the accuracy of the Cowling 
approximation and confirmed certain amount of errors of 
eigenfrequencies between them. Their results for rapidly ro- 
tating stars (Figure 4 in Jones et al. 2002) is well-matched 
with our result for rigid rotation case (Figure 1), although 
they have computed upper branch of f-modo and our re- 
sults are for lower branch. This similarity shows our result 
is consistent. 

There are two reasons why the accuracy of the Cowl- 
ing approximation becomes worse for differential rotation. 
The first reason is related to the density distribution of the 
equilibrium stars. In Figs. 3 and 4, equi-density contours of 
the rigidly and differentially rotating stars with N = 1.0 
and Tp/re = 0.6 are shown, respectively. Here rp and rc arc 
the polar radius and the equatorial radius, respectively. By 
comparing Fig. 3 and Fig. 4, in the rapidly and differentially 
rotating star, the dense region near the center extends to the 
outer region. Remember that the amplitude of the f-mode 
oscillation is significant in the outer region of the star. Then, 
in the differentially rotating star, the gravitational potential 
perturbation comes to play a relatively important role, be- 
cause the region with high-amplitude is dense. Therefore, 
the Cowling approximation which neglects the perturbed 
gravitational potential from the basic equations gives less 
accurate solutions. 

The second factor is the behavior of the eigenfunction 
due to differential rotation. In Figs. 5 and 6, the density per- 



turbation 5p is shown as a function of the relative distance 
from the rotation axis for rigidly (Fig. 5) and differentially 
(Fig. 6) rotating N = 1.0 polytropes. In the upper pan- 
els of these figures denote the eigenfunctions for the almost 
spherical configurations whose axis ratio is r-p/rc = 0.99, 
while the lower panels show those for the rapidly rotating 
stars with r-p/re = 0.60. The eigenfunctions in those figures 
are normalized by the maximum amplitude of the density 
perturbation in the star. From Fig. 5, the eigenfunctions of 
rigidly rotating polytropes are the monotonically increasing 
functions of the distance from the rotation axis and their am- 
plitudes become the maximum at the surface. This behaviors 
arc common with the previous results given in the cases of 
spherical stars and/or rigidly rotating stars. In Fig. 6, how- 
ever, the behavior of the eigenfunctions drastically change. 
When the stellar rotation is slow, the qualitative behaviors 
of the eigenfunctions that increase monotonically with the 
stellar radii are the same as those of rigidly rotating poly- 
tropes. When the stellar rotation is sufficiently rapid, the 
maximum values of the eigenfunctions appear in the inner 
region, not at the surface. This result means that in rapidly 
and differentially rotating stars the amplitudes of the modes 
in the high density region are rather large and that the ef- 
fects of the gravitational perturbation become important. 

In order to see the dependence of the accuracy of the 
Cowling approximation on the equations of state in differ- 
entially rotating stars, the eigenvalues of the m = 2 f-modes 
in differentially rotating {A = 1.0) stars, with several poly- 
tropic indices are also computed. The eigenvalues as func- 
tions of T/|VK| are shown in Fig. 7. Three curves of different 
types correspond to different models with N = 0.5, N = 1.0, 
and A'^ = 1.5. Similarly as shown in Fig. 1, the upper se- 
quences of the pair curves show the eigenvalues for the full 
perturbed systems, while the lower sequences correspond to 
the results of the Cowling approximation. 

In Fig. 8, the differences of the eigenfrequencies between 
the solutions of the full systems and those of the Cowling 
approximation are shown as in Fig. 2. For rigidly rotating 
stars, the high compressibility makes the accuracy of the 
approximation better, because the density distributions in 
highly compressible star is considerably condensed to the 
central region, so the density will be significantly low in the 
outer region where the amplitudes of the oscillations become 
laxge. In Fig. 7 and 8, almost the same tendency with the 
previous one can be observed. For N = 1.5 polytrope, the 
matter of the star is more centrally condensed compared 
with the smaller values of N. Hence, even if the star is ro- 
tating differentially with certain degree, the accuracy of the 
Cowling approximation does not become inaccurate. On the 
other hand, for N = 0.5 model, since the density of the star 
does not drop significantly in the outer region, the accuracy 
of the Cowling approximation becomes considerably inaccu- 
rate mainly due to the effect of differential rotation. 



4 DISCUSSION 
4.1 f-mode 

In the previous section, we have compared the results given 
by the system of the full perturbed equations with that in 
the Cowling approximations, and we have found that the 
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Figure 1. The eigenvalues of tiie m = 2 f-mode oscillations 
for JV = 1.0 polytropes. The eigenvalue a is normalized as 
a = (T/\/4nGpc- Three curves of different types (solid, dashed and 
dotted) denote sequences with rigid rotation, wealcly differential 
rotation (A = 1.0), and strongly differential rotation {A = 0.5), 
respectively. Tlic upper sequence of the curves of the same type 
(bold) corresponds to the solutions of the full perturbed equations 
and the lower sequence corresponds to those of the Newtonian 
Cowling approximation. 
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Figure 3. Density contours of a rigidly rotating star with the 
polytropic index N = 1.0 and rp/re = 0.6 in the meridional 
plane. The ith contour from the center denotes the isodensity 
curve with p = (l — O.lx i)pmax. Here pmax is the maximum 
density. 
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Figure 2. The differences of the eigenfrequcncies between the re- 
sults of the full perturbation and those of the Newtonian Cowling 
Approximation (5ct = Icrfun-pcrturb. — CapproxJ) for rigid rotation 
(solid), 1.0 (dashed), and A = 0.5 (dotted). The terminal point 
of the rigid rotation sequence corresponds to the mass shedding 
limit. 



accuracy of the approximation becomes worse when we con- 
sider rapidly and differentially rotating stars. And the de- 
gree of declination of the accuracy depends on the strength 
of differential rotation. Two reasons for this tendency can 
be considered. The first reason will be the prolongational 
effect of the equilibrium density distribution along equato- 
rial plane, by the presence of differential rotation. The sec- 
ond one will be the change that the maximum points of the 
eigenfunctions appear in the interior regions of the star, not 
surface regions. 

In this section, we will quantitatively estimate these ef- 
fects of differential rotation. For this purpose, we will com- 
pute the gravity perturbations numerically in Newtonian 
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Figure 4. Same as Fig. 3 but for differential rotation (A = 0.5). 

gravity. As mentioned above, since the reliability of the 

Cowling approximation becomes bettor when the density of 
the fluid components is sufficiently low in the region where 
the amplitudes of the oscillation get relatively larger, in or- 
der to estimate the relation between the density in an equi- 
librium and the amplitudes of the potential perturbations , 
we define a new quantity as follows: 



pS4>rn ex.-p{i{crt — rmp)) dV, 



(4) 



where, dtp = Sip/Scpmax is the perturbed gravitational poten- 
tial normalized by the maximum amplitude of 5<j) in the star. 

* These quantities, W^, for several sequences N = 1.0 
polytropes are shown in Fig. 9, for A = 1.0, and 0.5. The 

* Additionally, here, since it is assumed that the equilibrium con- 
figurations are axisymmetric and the perturbations have the de- 
pendency as exp(i(iTt — rmp)) in the aaimuthal direction, the in- 
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Figure 5. The eigenfunctions of the density perturbation for the 
m = 2 f-mode oscillations in the rigidly rotating polytropes with 
N = 1.0 and r-pjre = 0.99 (upper) and rp/re = 0.60 (lower). 
The horizontal axis denotes the relative distance from the rota- 
tion axis, r/rs{0), and the vertical axis is the normalized ampli- 
tude of the perturbed density. Here rs{0) denotes the surface of 
the unperturbed equilibrium configuration. Three different curves 
correspond to the eigenfunctions on the spokes oi = 7r/6,7r/3 
and tt/2 directions. 
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Figure 6. Same as Fig. 5 but for differential rotation (A = 0.5). 
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Figure 7. Same as Fig.l but for sequences with A = 1.0. Three 
curves of different types (solid, dashed and dotted) denote se- 
quences with N = 0.5, 1.0, and 1.5, respectively. 
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Figure 8. Same as Fig. 2 but for sequences with A = 1.0. Three 
curves of different types represent sequences with A'' = 0.5 (solid), 
N = 1.0 (dashed), and N = 1.5 (dotted). 



horizontal axis denotes r/|VK|. Here, actually plotted val- 
ues, W^, are normalized by the masses of the equilibrium 
models: 



Mo = j p 



dV. 



(5) 



Introducing such a quantity and comparing Figs. 9 with 
Fig. 2, we can see that the behavior of becomes consis- 
tent with the absolute values of 5a in terms of the function of 
A and T/|W^|. These tendencies of VK^ strongly suggest that 
our conjecture described in previous sections axe correct; i.e. 
the accuracy of the Cowling approximation is related with 
the values of W4,. 



tegral in the (^-direction, J exp(i[at — rmp)) dtp, need not to be 
carried out in the present context. In this meaning, we express 
this quantity as W^. 
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Figure 9. The values of for N = 1.0 polytropcs. Results for 
rigid rotation, weakly differential rotation {A = 1.0), and strongly 
diflferential rotation {A = 0.5) are shown by solid line, dashed line, 
and dotted line, respectively. 

4.2 r-mode 

In these days, from astrophysical point of view, r-mode os- 
cillations in neutron stars are noticed. A lot of studies which 
suggest that r-modes can be unstable in every rotating star 
and the accompanied gravitational waves can be detected by 
improved interferometer as like LIGO II has been done (An- 
dersson 1998, Friedman and Morsink 1998, for recent works, 
see the review of Andersson 2002 and its references). 

Since the eigenfunctions of the r-mode oscillations axe 
surface-standing, they may not suffer from the effect of ex- 
tension of the central dense regions. Additionally the gravity 
perturbation is much smaller than the fluid perturbations 
for r-modes (about the feature of r-mode oscillations, Saio 
1982, Yoshida and Lee 2002, for the eigenvalue of r-modes in 
differentially rotating stars, see Karino et al. 2001). Because 
of these reasons, the Cowling approximation may be accu- 
rate in the calculations of r-mode oscillations. In order to 
prove this, we have carried out the same computation with 
f-modcs for the m = 2 r-mode. The results are shown in Fig- 
ure 10. The plot is in the same manner with Figure 1. From 
this result, it can be recognized that the approximation will 
give good results for r-mode oscillations. This conclusion is 
also comparable with the assertion of Jones et al. (2002). 



5 CONCLUDING REMARKS 

In the present study, the accuracy of the Cowling approx- 
imation has been investigated for rapidly and differentially 
rotating stars in Newtonian gravity. We have adopted a poly- 
tropic equation of state with polytropic indices N = 0.5, 1.0, 
and 1.5 as the stellar models and estimated the accuracy of 
the approximation for the m = 2 f- and r-mode oscillations. 

For the f-mode in rigidly rotating stars, the tendency 
that the accuracy gets better for rapidly rotating stars has 
been reproduced. On the other hand, when the angular mo- 
mentum distribution of the star is not uniform, the accuracy 
of the Cowling approximation becomes worse when the ro- 
tation rate gets higher. The accuracy of the approximation 
also depends on the equation of state. For softer equations of 
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Figure 10. The same as Fig. 1. But in this case, the eigenvalues 
are divided by the central angular velocity Oc in order to compare 
the results with Karino et al. (2001). The lower curves in each pair 
denote the results obtained by the approximation. 

state, the accuracy of the approximation becomes relatively 
bc'ttor oveu if the effect of differential rotation exists, be- 
cause of the central condensation of the matter. Contrary, for 
stiffer equations of state, which are usually used to describe 
the neutron stellar matter, the accuracy of the approxima- 
tion becomes worse for rapidly and differentially rotating 
stars. The reasons why the accuracy becomes worse can be 
considered as follows. First, when the angular momentum 
distribution is not uniform, the equilibrium configurations 
arc highly deformed and the high density region extends to 
the outer part of the equatorial region. Second, when the 
angular momentum distribution is not uniform, the maxi- 
mum amplitude of the eigenfunction appears well inside the 
star where the equilibrium density is rather high. It is pos- 
sible to quantitatively confirm that the oscillations become 
large in the high density region, by using a new quantity, 
W4,, which is the integration of the density weighted by the 
eigenfunctions of the oscillations. 

As discussed in Introduction, it is extremely difficult 
to solve oscillations of general relativistic rotating compact 
stars. This is mainly because it is difficult to treat the per- 
turbations of the metric, i.e. the perturbed gravitational po- 
tentials. Thus it is desirable to apply the Cowling approxi- 
mation to general relativistic configurations and obtain rea- 
sonably accurate approximated solutions. 

Can we say anything about the applicability of the 
Cowling approximation to general relativistic configurations 
? It is widely believed that the accuracy of the Cowling ap- 
proximation is better for general relativistic configurations, 
because the strong gravity makes the effective central con- 
densation of the stellar matter higher. This effect of gravi- 
tation is similar to the effect of the increeising of the com- 
pressibility. As shown in the present paper, the ascent of 
the compressibility achieves good accuracy of the approxi- 
mation. Therefore, in GR, it is expected that the accuracy of 
the Cowling approximation would become better than that 
for Newtonian gravity. 

However, also in the case of GR, the same problem ac- 
companied with differential rotation will happen. When we 
apply the Cowling approximation in differentially rotating 
relativistic stars, the accuracy of the approximation will de- 
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pend on the positive effect of strong gravity and the nega- 
tive effect of differential rotation. The detailed study must 
be carried out to discuss about the accuracy of such an 
approximation in differentially rotating compact stars. As 
mentioned in the last section, however, the Cowling approx- 
imation works very well in the calculations of r-mode os- 
cillations (also see, Yoshida & Lee 2002). Hence, it is good 
idea to apply the approximation to obtain the eigenfrequen- 
cies of r-modes in general relativity. Also it is shown that 
the accuracy of the Cowling approximation becomes bet- 
ter for higher modes (Jackson 1995), therefore to compute 
those higher modes by using this approximation is one of 
the realistic usages of the Cowling approximation. 
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